AUBRY-MATHER MEASURES IN THE NON CONVEX SETTING* 

F. CAGNETTit, D. GOMESt, AND H. V. TRAN* 

Abstract. The adjoint method, introduced in [EvalO] and [Tra], is used to construct analogs 

to the Aubry-Mather measures for non convex Hamiltonians. More precisely, a general construction 

of probability measures, that in the convex setting agree with Mather measures, is provided. These 

measures may fail to be invariant under the Hamiltonian flow and a dissipation arises, which is 

■ ■ ■ described by a positive semi-definite matrix of Borel measures. However, in the case of uniformly 

' ^ , quasiconvex Hamiltonians the dissipation vanishes, and as a consequence the invariance is guaranteed. 
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1. Introduction. Let us consider a periodic Hamiltonian system whose energy 
^N ■ is described by a smooth Hamihonian iJ : T" x M" — > R. Here T" denotes the n- 

dimensional torus, n S N. It is well known that the time evolution 1 1— > (x(t), p(t)) of 
PLh ' the system is obtained by solving the Hamilton's ODE 



< 



±= -DpH{x,p), 

^ ■ < (1-1) 

p^ DxH{x,p). 



Assume now that, for each P G M", there exists a constant H{P) and a periodic 
function u{-,P) solving the following time independent Hamilton- Jacobi equation 

^ . H{x, P + D,u{x, P)) = H{P). (1.2) 

t^^ ■ _ 

Suppose, in addition, that both u{x, P) and H{P) are smooth functions. Then, if the 

ff^ . following relations 

X = x + Dpu{x,P), p = P + D,u{x,P), (1.3) 

^^ , define a smooth change of coordinates X{x^p) and P{x,p), the ODE (1.1) can be 

^— ^ ■ rewritten as 

(x=-DpH{P), 

X ^^ 

H , 

5t 1 Since the solution of (1.4) is easily obtained, solving (1.1) is reduced to inverting the 

change of coordinates (1.3). Unfortunately, several difficulties arise. 

Firstly, it is well known that the solutions of the nonlinear PDE (1.2) are not 
smooth in the general case. For the convenience of the reader, we recall the definition 
of viscosity solution. 
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Definition 1.1. We say that u is a viscosity solution of (1.2) if for each v € 
C°°(R") 

• If u — V has a local maximum at a point xq G M" then 

H{xo,P + Dvixo))<HiP); 

• If u — V has a local minimum at a point xq G R" then 

H{xo,P + Dvixo))>HiP). 

One can anyway solve (1.2) in this weaker sense, as made precise by the following 
theorem, due to Lions, Papanicolaou and Varadhan. 

Theorem 1.2 (See [LPV88]). Let iJ : T" x R" ^ M &e sm,ooth such that 

lim II{x,p) = +00. (1-5) 

|p|— S- + 00 

Then, for every P G M" there exists a unique II{P) G M such that (1.2) admits a Z"- 
periodic viscosity solution u(-,P) : T" ^ R. We call (1.2) the cell problem. It can 
be proved that all the viscosity solutions of the cell problem are Lipschitz continuous, 
with Lipschitz constants uniformly bounded in P. 

A second important issue is that the solution u(-, P) of (1.2) may not be unique, 
even modulo addition of constants. Indeed, a simple example is given by the Hamil- 
tonian H{x,p) — p ■ {p — Dip{x)), where ■)/' : T" — )■ M is a smooth fixed function. In 
this case, for P = and if (0) — 0, the cell problem is 

Du- D{u-il)) ==0, 

which admits both u = and u = ijj as solutions. Therefore, smoothness of u{x, P) 
in P cannot be guaranteed. 

Finally, even in the particular case in which both u{x,P) and II{P) are smooth, 
relations (1.3) may not be invertible, or the functions X{x,p) and P{x,p) may not be 
smooth or globally defined. 

Therefore, in order to understand the solutions of Hamilton's ODE (1.1) in the 
general case, it is very important to exploit the functions -ff(-P) and u{x,P), and to 
extract any possible information "encoded" in II{P) about the dynamics. 

1.1. Classical Results: the convex case. Classically, the additional hypothe- 
ses required in literature on the Hamiltonian H are: 
(i) -ff (x, •) is strictly convex; 

(ii) II{x,-) is superlinear, i.e. 

,. H{x,p) 
lim — -- — = +00. 
IpK+oo \p\ 

A typical example is the mechanical Hamiltonian 

H{x,p)^^^ + V{x), 

where V^ is a given smooth Z"-periodic function. Also, one restricts the attention to a 
particular class of trajectories of (1.1), the so-called one sided absolute minimizers of 
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the action integral. More precisely, one first defines the Lagrangian L : T" x M" — > M 
associated to H as the Legendre transform of H: 

L{x, v) := H*{x, v) = sup {-p ■ v - H{x,p)} for every (x, v) G T" x M". (1.6) 

Here the signs arc set following the Optimal Control convention (see [FS93]). Then, 
one looks for a Lipschitz curve x(-) which minimizes the action integral, i.e. such that 

i(x(i),x(i))dt< / L{y{t),y{t))dt (1.7) 



for each time T > and each Lipschitz curve y(-) with y(0) ~ x(0) and y(T) = x(T). 
Under fairly general conditions such minimizers exist, are smooth, and satisfy the 
Euler-Lagrange equations 

^ [D,L{^{t),±{t))] - Z?,i(x(i), i(t)), t e (0, +(»). (1.8) 

It may be shown that if x(-) solves (1.7) (and in turn (1.8)), then (x(-),p(-)) is a 
solution of (1.1), where p(-) := — I?^L(x(-),x(-)). This is a consequence of assump- 
tions (i) and (ii), that in particular guarantee a one to one correspondence between 
Hamiltonian space and Lagrangian space coordinates, through the one to one map 

$ : T" X M" ^ T" X M" defined as 

$(x, v) :== (x, -D^L{x, v)). (1.9) 

There are several natural questions related to the trajectories x(-) satisfying (1.7), 
in particular in what concerns ergodic averages, asymptotic behavior and so on. To 
address such questions it is common to consider the following related problem. 

In 1991 John N. Mather (see [Mat91]) proposed a relaxed version of (1.7), by 
considering 



mm 



L{x,v) di>{x,v), (1-10) 



T"xl 



where V is the class of probability measures in T" x R" that are invariant under 
the Euler-Lagrange flow. In Hamiltonian coordinates the property of invariance for a 
measure v can be written more conveniently as: 

{(t),H}dfi{x,p) = 0, for every e C](T" x M"), 

T"xK" 

where /i = $#^ is the push-forward of the measure v with respect to the map $, i.e., 
the measure fx such that 



(l){x,p) dii{x,p) — / (j){x,—DyL{x,v))dh'{x,v), 

T"xE" JT"xR" 

for every (j) E Cc(T" x M"). Here the symbol {•, •} stands for the Poisson bracket, that 
is 

{F, G} :^ DpF ■ D^G - D^F ■ DpG, for every F,G e G\T' x M"). 
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Denoting by 7'(T" x M") the class of probability measures on T" x M", we have 
V= liye T'CT" X M") : [ {(j), H} d<^^v{x,p) ^ 0, for every (f) € C,^(T" x M") I 

(1.11) 

The main disadvantage of problem (1.10) is that the set (1.11) where the mini- 
mization takes place depends on the Hamiltonian H and thus, in turn, on the integrand 
L. For this reason, Ricardo Mane (see [Mn96]) considered the problem 



min / L{x,v) di^{x,v), (1-12) 

i^GJ^Jt"xR" 



where 



J" := <iye 7'(T" x R") : / v D%l){x) dv{x,v) = 0, for every V- e C^(T") 

I JT"xR" 

Measures belonging to T are called holonomic measures. Notice that, in particular, 
to every trajectory y(-) of the original problem (1.7) we can associate a measure 
^'y(.) G T. Indeed, for every T > we can first define a measure ^'T.y() G V{T" x M") 
by the relation 

1 '■^ 



x,v)duT.y{-){x,v):^- <j){y{t),y(t))dt for every G Cc(T" x M"). 

T"xR" ^ JO 

Then, from the fact that 

snppi^T.yi) C T"x[-M, M], for every T > 0, (M = Lipschitz constant of y(-)) 

we infer that there exists a sequence Tj — )■ oo and a measure Uyf.^ E 7'(T" x R") such 
that VT-,y(-) ~^ ^y{-) in^ the sense of measures, that is, 

1 '•^' 



lim — / (j){y{t),y{t))dt^ (j){x,v) dvyi.){x,v) for every </> G Cc(T"xR" ). 

j^oo Tj Jq Jt"xR" 

(1.13) 
Choosing (l){x,v) — v ■ Dtp{x) in (1.13) it follows that fy(.) G T, since 

vDi;{x)diyy^.)ix,v)^ lim 1 / ' yityD^P{yit)) dt ^ lim tkiM_^tkM^o, 

In principle, since !F is much larger than the class of measures T), we could expect 
the last problem not to have the same solution of (1.10). However, Mane proved that 
every solution of (1.12) is also a minimizcr of (1.10). 

A more general version of (1.12) consists in studying, for each P G R" fixed, 

min / {L{x^v) -'r P ■ v) dv{x,v), (1-14) 

''6-^JT"xR" 

referred to as Mather problem. Any minimizer of (1.14) is said to be a Mather measure. 
An interesting connection between the Mather problem and the time independent 
Hamilton- Jacobi equation (1.2) is established by the identity: 

-H{P)^ mm {L{x,v) + P-v) diy{x,v). (1.15) 



Notice that problems (1.12) and (1.14) have the same Euler-Lagrange equation, but 
possibly different minimizers, since the term P ■ v is a null Lagrangian. The following 
theorem gives a characterization of Mather measures in the convex case. 

Theorem 1.3. Let iJ : T" x R" — >• R &e a smooth function satisfying (i) and 
(a), and let P G E". Then, v G ^(T" x M") is a solution of (1.14) if and only if: 

(a) / H{x,p)dii{x,p) = H{P) = H{x,p) ii-a.e.; 

JT"xK" 

(6) / {p-P)-DpH{x,p)d^{x,p)^0; 

JT"xK" 

(c) f DpH{x,p)-D(f){x)dn{x,p)^0, for every <j) e C\T") , 

JT"xE" 

where /i ~ 'I'#^ (md H{P) is defined by Theorem 1.2. Before proving Theorem 1.3 
we state the following proposition, which is a consequence of the results in [Mn96], 
[Fat97a], [Fat97b], [Fat98a], [Fat98b] and [EGOl]. 

Proposition 1.4. Let H : T" x R" -> M 6e a smooth function satisfying (i) and 
(a). Let P e R", let v G 7^(1" x M") he a minimizer of (1.14) and set fi = $#i^. 
Then, 

(1) n is invariant under the Hamiltonian dynamics, i.e. 

{</), H} d^i{x, p)^0 for every (j) € C\ (T" x M"); 

T"xK" 

(2) /i is supported on the graph 

S :={(x,p) e T" X M" -.p^P + D^uix^P)}, 

where u is any viscosity solution of (1.2). We observe that property (2), also known 
as the graph theorem, is a highly nontrivial result. Indeed, by using hypothesis (ii) 
one can show that any solution u{-,P) of (1.2) is Lipschitz continuous, but higher 
regularity cannot be expected in the general case. 

Proof. [Proof of Theorem 1.3] To simplify, we will assume P = 0. 

Let i^ be a minimizer of (1.14). By the previous proposition, we know that 
properties (1) and (2) hold; let us prove that /i = ^#i^ satisfies {a)~{c). By (1.15), 
we have 

/ L{x, v) dv{x, v) = -ff (0). 

JT"xR" 

Furthermore, because of (2) 

H{x,p)d^l{x,p)^H{{)), 

/T"xR" 

that is, (a). Since H{x,p) = —L{x,—DpH{x,p)) +p- DpH{x,p), this implies that 

p ■ DpH{x,p) d^i{x,p) = 0, 

T"xR" 

and so (b) holds. Finally, (c) follows directly from the fact that v ^ T . 
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Let now /i e 7'(T" x K") satisfy {a)-{c), and let us show that u ~ {^ ^)#l^ is a 
minimizer of (1.14). First of all, observe that v ^ T . Indeed, by using (c) for every 

^ g Ci(T") 



V ■ D')p{x) dv{x, v) = — DpH{x,p) ■ Dtplx) dji{x,p) — 0. 

T"xR" JT"xR" 

Let now prove that :^ is a minimizer. 

Integrating equality H{x,p) ~ —L{x,—DpH{x,p)) + p ■ DpH{x,p) with respect 
to fx, and using (a) and (b) we have 



HiO)= f H{x,p)dii{x,p) 

JT"xR" 



L{x,-DpH{x,p))d^{x,p) + / p- DpH{x,p)d^{x,p) 

L{x,—DpH{x,p))dij{x,p) = — / L{x,v) di^{x,v). 

T"xR" JT"xR" 

By (1.15), z^ is a minimizer of (1.14). 
D 

1.2. The Non Convex Case. The main goal of this paper is to use the tech- 
niques of [EvalO] and [Tra] to construct Mather measures under fairly general hy- 
potheses, when the variational approach just described cannot be used. Indeed, when 
(i) and (ii) are satisfied H coincides with the Legendre transform of L, that is, identity 
H = H** holds. Moreover, L turns out to be convex and superlinear as well, and 
relation (1.9) defines a smooth diffeomorphism, that allows to pass from Hamiltonian 
to Lagrangian coordinates. 

First of all, we extend the definition of Mather measure to the non convex setting, 
without making use of the Lagrangian formulation. 

Definition 1.5. We say that a measure ^ £ V{T" x M") is a Mather measure if 
there exists P G R" such that properties (a)-(c) are satisfied. The results exposed in 
the previous subsection show that, modulo the push-forward operation, this definition 
is equivalent to the usual one in literature (see e.g. [Fat], [Mn96], [Mat91]). We would 
like now to answer the following natural questions: 

• Question 1: Does a Mather measure exist? 

• Question 2: Let /x be a Mather measure. Are properties (1) and (2) satisfied? 
We just showed that in the convex setting both questions have affirmative answers. 
Before addressing these issues, let us make some hypotheses on the Hamiltonian H. 
We remark that without any coercivity assumption (i.e. without any condition similar 
to (ii)), there are no a priori bounds for the modulus of continuity of periodic solutions 
of (1.2). Indeed, for n = 2 consider the Hamiltonian 

H{x,p) ~ p\— p\ for every p — {pi,p2) G R^. 

In this case, equation (1.2) for P = and H{P) = becomes 

ul-ul^O. (1.16) 

Then, for every choice of / : R ^ R of class C^, the function u{x,y) = f{x — y) is 
a solution of (1.16). Clearly, there are no uniform Lipschitz bounds for the family of 
all such functions u. Throughout all the paper, we will assume that 
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(HI) H is smooth; 

(H2) H{-,p) is Z"-pcriodic for every p € M"; 

(H3) lini|p|^_|_oo — 1-^(2;, p)P + DxH{x,p) ■ p] = +00 uniformly in x. 

Note that if hypothesis (ii) of the previous subsection holds uniformly in x and we 
have a bound on DxH{x,p), e.g. \DxH{x,p)\ < C(l + |p|), then (H3) holds. 

First we consider, for every £ > 0, a regularized version of (1.2), showing existence 
and uniqueness of a constant H (P) such that 



-Au^ix) + H{x, P + Du^{x)) ^ H'(P) 



(1.17) 



admits a Z"-periodic viscosity (in fact smooth) solution (see Theorem 2.1). 

Thanks to (H3), we can establish a uniform bound on ||I?w^||loo and prove that, 
up to subsequences, H (P) -^ H{P) and u^{-,P) converges uniformly to u(-,P) as 
e -^ 0, where H{P) and u{-,P) solve equation (1.2). 

Observe that, in particular, this shows that Theorem 1.2 still holds true under 
assumption (H3) when (1.5) does not hold, as for instance when n = 1 and 

H{x,p) = p^ + V{x)^ V smooth and Z"-periodic. 

On the other hand (1.5) does not imply (H3), sec the Hamiltonian 

H{x,p) =p^ f3 + sin(eP''(cos27ra;))' 



(here again n = 1). Thus, although (H3) seems to be a technical assumption strictly 
related to the particular choice of the approximating equations (1.17), it is not less 
general than (1.5), as just clarified by the previous examples. Anyway, it is not clear at 
the moment if the results we prove in the present paper are still true for Hamiltonians 
satisfying (1.5) but not (H3). 

Once suitable properties for the sequence {w^} are proved, for every e > we 
define the perturbed Hamilton SDE (see Section 3) as 



dp" = D^H{yf ,p") dt + eD'^u" dwt, 



(1.18) 



where Wt is a n-dimensional Brownian motion. The main reason why we use a stochas- 
tic approach, is that in this way we emphasize the connection with the convex setting 
by averaging functions along trajectories. Nevertheless, our techniques can also be 
introduced in a purely PDE way (see Section 3.3 for a sketch of this approach). 

In the second step, as just explained, in analogy to what is done in the convex setting 
we encode the long-time behavior of the solutions t h^ (x^(t),p^(t)) of (1.18) into a 
family of probability measures {/i'^}e>Oi defined by 



(f){x,p) d^"{x,p) :— lim 



T"xl 



1 

t' 



■E 



^{x'{t),p'{t))dt 



for every (^e Cc(T"xIR"), 



where with E[-] we denote the expected value and the limit is taken along appropriate 
subsequences {Tj}jgN (see Section 3.1). 

Using the techniques developed in [EvalO], we are able to provide some bounds on 
the derivatives of the functions u^. More precisely, defining 0^<s as the projection on 
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the torus T" of the measure ^f (see Section 3.2), we give estimates on the (L^, (i0^»)- 
norm of the second and third derivatives of u'^, uniformly w.r.t. e (see Proposition 
4.1). 

In this way, we show that there exist a Mather measure ^ and a nonnegative, 

symmetric n x n matrix of Borel measures {mkj)k.j=i n such that fj,'^ converges 

weakly to fi up to subsequences and 

{0, H}d^i+ f <j)p^p^ dnikj =0, V</> G C2(T" X R"), (1.19) 

T"xE" JT"xR" 

with sum understood over repeated indices (see Theorem 5.1). As in [EvalO], we call 
nikj the dissipation measures. Relation (1.19) is the key point of our work, since it 
immediately shows the differences with the convex case. Indeed, the Mather measure 
fi is invariant under the Hamiltonian flow if and only the dissipation measures rukj 
vanish. When H{x,-) is convex, this is guaranteed by an improved version of the 
estimates on the second derivatives of u'^ (see Proposition 4.1, estimate (4.4)). We 
give in Section 10 a one dimensional example showing that, in general, the dissipation 
measures {mkj)k,j=i,....n do not disappear. 

We study property (2) in Section 8. In particular, we show that if (1.2) ad- 
mits a solution u{-,P) of class C^, which is a rather restrictive condition, then the 
corresponding Mather measure /i given by Theorem 1.2 satisfies 

DpH{x, P + D,u{x, P))-{p-P- D,u{x, P)) = 

in the support of /i (see Corollary 8.2). Observe that this single relation is not enough 
to give us (2) in general, e.g. n > 2. 

Finally, we are able to provide some examples of non-convex Hamiltonians (see 
Section 9), for which both properties (1) and (2) are satisfied. We observe that the 
case of strictly quasiconvex Hamiltonians, which appears among our examples, could 
also be studied using duality (see Section 9.7). 

2. Elliptic regularization of the cell problem. We start by quoting a clas- 
sical result concerning an elliptic regularization of equation (1.2). This, also called 
vanishing viscosity method, is a well known tool to study viscosity solutions. In the 
context of Mather measures this procedure was introduced by Gomes in [Gom02] , see 
also [Ana04], [AIPSM05], [ISM05]. 

Theorem 2.1. For every e > and every P G M", there exists a unique number 
H (P) G K such that the equation 

l^u^{x) + H{x, P + Du^ix)) = H\P) (2.1) 

admits a unique (up to constants) 17^ -periodic viscosity solution. Moreover, for every 
PG M" 

lim H (P) = H{P), and u^ ^>- u uniformly (up to subsequences), 

where H{P) G M and u : T" — > R are such that (1.2) is satisfied in the viscosity sense. 
We call (2.1) the stochastic cell problem. 

Definition 2.2. Lete > andP G W\ The linearized operator L'^-^ : C^(T") -^ 
C(T") associated to equation (2.1) is defined as 

p2 

L^^^v{x) := Av{x) + DpH{x, P + Du%x)) ■ Dv{x), 



for every v E C^(T"). 

Proof. [Sketch of the Proof] We mimic the proof in [LPV88]. For every A > 0, 
let's consider the following problem 

Xv^ + H{x, P + Dv^) = —Av^. 

The above equation has a unique smooth solution v^ in M" which is ^"-periodic. 
We will prove that || Au'^H^oo , ||Dw^||l=° < C, for some positive constant C independent 
on A and e. By using the viscosity property with (p = as a test function, we get 

1 1 Aw^ 1 1 Loo < C. Let now w^ = . Then we have 

2 2 

2\w^ + DpH ■ Dw^ + D^H ■ Dv^ = —Aw^ - —\D^v^\^. 
Notice that for e < 1 





2 ' ' - 4 ' ' ^ ' - 2 


Therefore, 






1 p2 

2\w^ + DpH ■ Dw^ + D^H ■ Dv^ + -H^ -C< y Aw^ 



At Xi e T" where 'w^{xi) — maxx" w^ 

2Xw^{xi) + D^H ■ Dv^(xi) + -H^ < C. 

Since w^(xi) > 0, using condition (H3) we deduce that w^ is bounded independently 

of A, e. Finally, considering the limit A — > we conclude the proof. 

D 

Remark 2.3. Bernstein method and (H3) were used in the proof to deduce the 
uniform hound on \\Dv^\\[j<x, which is one of the key properties we need along our 
derivation. See [Lio82, Appendix 1] for conditions similar to (H3). 

The classical theory (see [Lio82]) ensures that the functions u^{-,P) are C°°. In 
addition, the previous proof shows that they are Lipschitz, with Lipschitz constant 
independent of e. 

3. Stochastic dynamics. We now introduce a stochastic dynamics associated 
with the stochastic cell problem (2.1). This will be a perturbation to the Hamilto- 
nian dynamics (1.1), which describes the trajectory in the phase space of a classical 
mechanical system. 

Let (T",(T, P)be a probability space, and let wt be a n-dimensional Brownian 
motion on T". Let e > 0, and let u^ be a Z"-periodic solution of (2.1). To simplify, 
we set P = 0. Consider now the solution x'^(i) of 

{dyf = -DpH{yf,Du^{yi^))dt + edwt, 

\x-(0)=x, ^-^ 

with X G T" arbitrary. Accordingly, the momentum variable is defined as 

l>%t)^Du^{^-{t)). 
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Remark 3.1. From Remark 2.3 it follows that 



sup|p'^(i)| < oo. 
t>o 

Let us now recall some basic fact of stochastic calculus. Suppose z : [0, +oo) 
is a solution to the SDE: 



rfzj = Oidt + hij wj 



i-1, 



,n, 



with ai and bij bounded and progressively measurable processes. Let ip : R" x . 
be a smooth function. Then, ip{z,t) satisfies the Ito formula: 

1. 



dip = (^2, dzi + iipt + -^bijbjkifz.zk 1 dt. 
An integrated version of the Ito formula is the Dynkin's formula: 



(3.2) 



E[^{z{T))-^{z{0))]^E 



^a,Dz^^{z{t)) + h,,b,kDi^,J(z{t)) ] dt 



Here and always in the sequel, we use Einstein's convention for repeated indices in a 
sum. In the present situation, we have 

ai = -Dp^H{x.^,Du''), bij = e^^ . 

Hence, recalling (3.1) and (3.2) 

dpi = u%.^. dx'j + -—A{u%.)dt = —L'^'^u%.dt + eu%.^. dwl 

^ D^Mdt + eul^^^dwi, 

where in the last equality we used identity (4.9). Thus, (x^,p'^) satisfies the following 
stochastic version of the Hamiltonian dynamics (1.1): 



dx.'^ = —DpH{x^,p'^) dt + edwt, 
dp" ^ D^H{x.",p") dt + eD^u" dwt 



(3.3) 



We are now going to study the behavior of the solutions u" of equation (2.1) along 
the trajectory x'^(i). Thanks to the Ito formula and relations (3.3) and (2.1): 

^2 



du"{x"{t)) = Du'dx" + y Au^ dt = -L^'^u^dt + sDu" dwt 



H - H - Du" ■ DpH ] dt + eDu" dwt. 
Using Dynkin's formula in (3.4) we obtain 



(3.4) 



H - H - Du" ■ Dr,H dt 



E{u%x%T))-u'{x.''{0))) ^E 



We observe that in the convex case, since the Lagrangian L is related with the Hamil- 
tonian by the relation 



L=p- DpH - H, 



we have 



u'{x'{Q))^E 



(L + H )dt + u"{x.%T)) 
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3.1. Phase space measures. We will encode the asymptotic behaviour of the 
trajectories by considering ergodic averages. More precisely, we associate to every 
trajectory (x'^(-),p'^(-)) of (3.3) a probability measure ^^ E 7'(T" x R") defined by 



/ 



T"xR" 



{x,p) dfj,%x,p) -.^ \im —E 

T— !-oo 1 



<l^{^%t),p%t))dt 



(3.5) 



for every (j) <E Cc(T" x R"). In the expression above, the definition makes sense 
provided the limit is taken over an appropriate subsequence. Moreover, no uniqueness 
is asserted, since by choosing a different subsequence one can in principle obtain a 
different limit measure fi^. Then, using Dynkin's formula we have, for every € 

C2(T" xK"), 



i?[0(x-(r),p-(r))-0(x-(o),p^(o))] = i? 



E 



T .^2 



Dp(j) ■ D^H - D^(j) ■ DpH I dt 



>x,x, + e^u" 



XiXjVXiPj + 2 ^XiXk'^XiX-jVpkPj I "''' 



(3.6) 



Dividing last relation by T and passing to the limit as T ^- +oo (along a suitable 
subsequence) we obtain 






£2 ^2 

'77'PxiXi +e "UxiXjVxiPj + 'TT'^^XiXk'^XiXj 



"PkPj 



dfi"" = 0. 

(3.7) 

3.2. Projected measure. We define the projected measure 9^^ G 7'(T") in the 
following way: 



Lp{x) dO ^^E {x) 



(p{x)dfi%x,p), 



V(/9eC(T"). 



Using test functions that do not depend on the variable p in the previous definition 
we conclude from identity (3.7) that 



DpH ■ D^ de^. = — 



AifidOp 



V93€C-'(T"). (3.: 



T"xl 



3.3. PDE Approach. The measures ^^ and 9pE can be defined also by using 
standard PDE methods from (3.8). Indeed, given m^ we can consider the PDE 

£2 

— A0' + div {DpH{x, Du") 6^) = 0, 

which admits a unique non-negative solution 9^ with 

9^{x)dx = 1, 

since it is not hard to see that is the principal eigenvalue of the following elliptic 
operator in C'^iT'"-): 



-Av - div{DpH{x, Du^) v). 
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Then fj,'^ can be defined as a unique measure such that 

tp{x,p) dii^{x,p) — / tp{x,Du'^{x)) dd'^{x), 

T"xR" JT" 

for every t/j E Cc(T" x M"). Finally, identity (3.7) requires some work but can also be 
proved in a purely analytic way. 

4. Uniform estimates. In this section we derive several estimates that will be 
useful when passing to the limit as e ^- 0. We will use here the same techniques as 
in [EvalO] and [Tra]. 

Proposition 4.1. We have the following estimates: 









(4.1) 


e^ [ \Dl^u-\^ d0^. < [ 


\Dpu'\^de^.+ f \DpH- 


-DpH'\'^d0f,., 


(4.2) 


e' f \Dul^,fd9,.<C 




i = 1, . . . ,n. 


(4.3) 



In addition, if H is uniformly convex inp, inequalities (4.1) and (4.2) can be improved 
to: 

f \Dl,u-\' d9^. < C, (4.4) 

/ \DJ,^u'fdef,.<CtTiicc{DlpH'), (4.5) 

respectively. Here C denotes a positive constant independent of e. 

Remark 4.2. Estimate (4.4) was already proven in [EvalO] and [Tra]. 

To prove the proposition we first need an auxiliary lemma. In the following, we 
denote by j3 either a direction in M" (i.e. /3 € M" with |/3| = 1), or a parameter (e.g. 
(3 — Pi for some i G {1, . . . , n}). When (3 — Pi for some i G {1, . . . , n} the symbols 
Hfi and -ff/3/3 have to be understood as Hp. and -ffp.p. , respectively. 

Lemma 4.3. We have 

e^ f \D^ulfd9^^ ^2 f u^Hl-Hp)de^e, (4.6) 

/ (hI^ ~ Hpn ^ 2DpHf3 ■ D^v'p - DlpHD^ifp ■ D^u%) d6^. - 0, (4.7) 

£2 / iD^u^^p de^,. = 2 / u'^p{Hl^ - Hpp - 2DpHp ■ D^u'^ - D^pH : D^u^p ® D^u^p dO^e. 

(4.8) 

Proof. By differentiating equation (2.1) with respect to /? and recalling Definition 
2.2 we get 

L^-Pu-p^Hl-Hp, (4.9) 

so that 
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Integrating w.r.t. 6*^= and recalling (3.8) we get (4.6). 

To prove (4.7), we differentiate (4.9) w.r.t. f3 obtaining 

L^'^M^^ = H^^ - Hpp - 2DpH[3 ■ D^ifp - DI^H : D,u% ® D^u^. (4.10) 

Integrating w.r.t. 9^^ and recalling (3.8) equality (4.7) follows. Finally, using (4.10) 



-Z/^' (JM^^I )—u'j)pL^' u^^ ^l-CxU^/sl 



^/3/3 



e |2 



{Hpp ~ Hf3(3 " 2DpH(3 ■ D^u'is - D'^pH : D^u'^s (g) D^m^) - — |D^w^^| 



Once again, we integrate w.r.t. 6^e and use (3.8) to get (4.8). D We can now proceed 
to the proof of Proposition 4.1. 

Proof. [Proof of Proposition 4.1] 

Summing up the n identities obtained from (4.6) with /3 — xi, . . . ,Xn respectively, 
we have 






D^u" ■ D^Hde,, 



Thanks to Remark 2.3, (4.1) follows. Analogously, relation (4.2) is obtained by sum- 
ming up (4.6) with /3 = Pi, P2, . . . , P„, which yields 



e^ / iDj^^u"]^ dOf^e =2 / Dpu" 



DpH -Dr,H 



de^.. 



Let us show (4.3). Thanks to (4.^ 
\D.u'fd0^, 



= -2 / ul^^^iH^^^^ + 2DpH^^ ■ D^ul^ + D^pH : D^u^^ «) D^uljdd^ 



Since the functions m^ are uniformly Lipschitz, we have 



\HxiXi\, \DpHxi\, \DppH\ < C, on the support of 6*^^. 



Hence, 



\Dxui^,fde^.<c 



<C I 



\DlX\d9^^ 



T" 



I^L"1'rfV 



iDLu^i^de,, 



\DtX\U9^, 



Finally, assume that H is uniformly convex. Thanks to (4.7) for every i = 1, . . . , n 

= / {Hx^x, + 2DpHx^ ■ Dxul^ + DlpHDxul^ ■ Dxuljde^e 

> / {Hx^x, + 2DpHx, ■ Dxul^) dOp. + a\\Dxul^\\l^„„.a 



:de„e)' 
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for some a > 0. Thus, using Cauchy's and Young's inequalities, for every r; e R 

Finally, 

[a — i] )\\D3:ul.\\^2(Y^.dg^^-f < HxiXidd/i' H ^-ll-Dp-ffx. ||i2(Tn;de^^)• 

Choosing rj^ < a we get (4.4). 

Let i e {1, . . .,n} and let us integrate w.r.t. 0^= relation (4.10) with /3 ~ Pf. 



0= [ (Hp^p^ - Hp^p^ - 2DpHp^ ■ Dxu'p^ - DlpHD^u^ ■ D^u^Jde^e . 

Since Dp H is positive definite, 

a f \Dxu'pf de^. < f {Hp^p^ - Hp^p^ - 2DpHp^ ■ DxU%)d9^. 
< / (H%^p^-2DpHp^-Dxu'pJd9f,e. 

Using once again Cauchy's and Young's inequalities and summing up with respect to 
i = 1,. . .,n (4.5) follows. D 

5. Existence of Mather measures and dissipation measures. We now look 
at the asymptotic behavior of the measures ^^ as e — > 0, proving existence of Mather 
measures. The main result of the section is the following. 

Theorem 5.1. Let _ff : T" x R" -^ K. be a smooth function satisfying conditions 
(H1)-(H3), and let {/i'^}£>o he the family of measures defined in Section 3. Then there 
exist a Mather measure fi and a nonnegative, symmetric n x n matrix {mkj)k,j=i,...n 
of Borel measures such that 

^^ ^^ fj, in the sense of measures up to subsequences, (5-1) 

and 

f {c^,H}dti+ f (f>p,p^ drukj ^ 0, V<^eC2(T"xR"). (5.2) 

Moreover, 

supp/i and supp TO are compact. (5-3) 

We call the matrix nikj the dissipation measure. 

Proof. First of all, we notice that since we have a uniform (in e) Lipschitz estimate 
for the functions u*^, there exists a compact set K C T^ x M" such that 

supp^" d K, Ve > 0. 
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Moreover, up to subsequences, we have (5.1), that is 



hm 

e-s-O 



^d^^ 



T"xR" 



^d^i, 



T"xl 



for every function (f) E Cc(T" x K"), for some probabihty measure ^ G 7'(T" x R"), 
and this proves (5.1). From what we said, it foUows that 

supp ^ C K. 

To show (5.2), we need to pass to the hmit in relation (3.7). First, let us focus on the 
second term of the aforementioned formula: 



T^xR" 



r2 £2 

-Z-<PxiXi +£ '^XiXjVxiPj + '7r''^XiXk'^XiXjVpkPj 



dn' 



(5.4) 



By the bounds of the previous section. 



lim 



T"xl 






dfi" = 0. 



However, as in [EvalO], the last term in (5.4) does not vanish in the limit. In fact, 
through a subsequence, for every k, j = 1, . . . , n we have 



T^xR" 



ul,xk'^l,x,i^{x,P) dfi'ix^p) 



ip(x,p)dmkj{x,p) W e Cc(T"xR"), 



for some nonnegative, symmetric nxn matrix {mkj)k.j=i,...n of Borel measures. Pass- 
ing to the limit as e — > in (3.7) condition (5.2) follows. From Remark 3.1 we infer 
that suppTTi C K, so that (5.3) follows. 

Let us show that n satisfies conditions (a)-(c) with P = 0. As in [EvalO] and 
[Tra], consider 

(h{x,p)-H'Y df,%x,p)^'- [ \Au%x)\Ufi%x,p)^0 

T"xR" ^ ' 4 Jt"xR" 

as £ —7- 0, where we used (2.1) and (4.1). Therefore, (a) follows. Let us consider 
relation (3.7), and let us choose as test function = (p{u'^). We get 



= 0. 



Passing to the limit as e — > 0, we have 

/ ip'{u)p- DpHdfi^O. 

JT"xR" 

Choosing ip{u) = u we get (b). Finally, relation (c) follows by simply choosing in (5.2) 

test functions that do not depend on the variable p. 

D 

We conclude the section with a useful identity that will be used in Section 9. 

Proposition 5.2. For every A G M 



JT"xR" 



dm 



kj 



0. 



(5.5) 
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Proof. First recall that for any function / : R ^ R of class C^ 

{H,f{H)} = 0, 

and, furthermore, for any tp G C^(T" x M") 

{H,^Pf{H)}^{H,i:}f{H). 

Let now A € M. By choosing in (5.2) (j) — ipf{H) with f{z) — e^^ and ip = 1 vfe 
conclude the proof. D 

6. Support of the dissipation measures. We discuss now in a more detailed 
way the structure of supp m. 
Proposition 6.1. We have 



suppm CZ 11 coG(x) =: K, (6.1) 

where with coG{x) we denote the convex hull in M" of the set G{x), and 
G(x) ■.^suppnn{{x,p) e T" xM" :a;=:r}, x e T". 

Remark 6.2. We stress that the convex hull of the set G{x) is taken only with 
respect to the variable p, while the closure in the right-hand side of (6.1) is taken in 
allT' X M". 

Proof. [Sketch of the proof] For t > sufficiently small, we can choose an open set 
K^ in T" X M" such that K C K^, dist (dKr, K) < r, and Kr{x) := {p e R" : (x,p) e 
Kt} is convex for every x e T". 

Also, we can find a smooth open set K2t C T" x R" such that, for every x € 
T", K2t{x) :— {p (z R" : {x,p) G K2t} is strictly convex, K2t{x) D Kr{x), and 
dist(ai4:2r(a;),iirr(x)) < r. 

Finally, we can construct a smooth function r/r : T" x R" -^ R such that for every 
xeT": 

• rjr{x,p) — for p G Kr{x). 

• p H- > rj-rix^p) is convex. 

• p I— > rjr{x,p) is uniformly convex on R" \ K2t{x). 
In this way, r]r{x,p) — on K^- D K D supp/i. Therefore 

/ {rjr,H}dfi^O. 

Combining with (5.2), 

JT"xR" 

which implies supp?Ti C [J^^jn K2t{x). Letting r — > 0, we finally get the desired 
result, n As a consequence, we have the following corollary. 
Corollary 6.3. 



suppm C co{H(x,p) < H}. 

Proof The proof follows simply from the fact that for every x € T" we have 

Gix)c{H{x,p)<H}. 

D 
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7. Averaging. In this section we prove some additional estimates concerning 
averaging with respect to the process (1.17). When necessary, to avoid confusion we 
will explicitly write the dependence on P. Let us start with a definition. 

Definition 7.1. We define the rotation number po associated to the measures fx 
and m as 



Po '■= lim lini E 



x^(T) - x=(0) 



T 



where the limit is taken along the same subsequences as in (3.5) and (5.1). The 
following theorem gives a formula for the rotation number. 

Theorem 7.2. 

There holds 



PO 



T"xl[ 



DpH dp,. 



Moreover, defining for every £ > the variable X*^ :— x*^ + Dpu^{x.^), we have 



E 



T 



= -DpH^{P), 



(7.1) 



(7.2) 



lim E 

T^ + oo 



X'^(r) - X'^(O) + DpH {P)T 

f 



<2ne^ + 2 / {Dpu^^ d0f,e 



2 / \DpH- DpH^l^dOf,. 



Proof. Choosing (j}{x) — xi with i = 1, 2, 3 in (3.6) we obtain 



E 



yf{T)-yf{Q) 



T 



-E 



^J DpH{^%t),p%t))dt 



Passing to the limit as T — > +oo 



Pe 



lim E 



x'^(T) - x'^(O) 



T 



DpHdp^. 



T"xl 



We get (7.1) by letting e go to zero. 

To prove (7.2), recalling Ito's formula (3.2) we compute 



d^^^d.^ + Di.^u^.nd.^ + -DpAu^.ndt 



-DpHix.", p'){I + Dj.^u'ix.")) + —DpAu^ix.") J dt + e{I + Df^^u^x.")) dwt, 

where in the last equality we used (3.1). By differentiating equation (2.1) w.r.t. P we 
obtain 

^DpHi^^, p^)(/ + D%Xi^')) + '^DpAu^i^') = ^DpH\p), (7.3) 
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so that 



dX^ = -DpH (P) dt + e{I + Dj^^u^x^)) dw- 



(7.4) 



Using the fact that 



E 



e{I + Dj,X{yi''))dwt 



-0, 



(7.2) follows. 

Finally, using once again Ito's formula (3.2) and relation (7.4) we can write 



eM2 



X.%t)-X.^O) + DpH {P)t 
= 2 (X.%t) - X'^(O) + DpH'{P)t\ {dX^ + DpH\P) dt) + e^\I + Dl^u'{^ 
= 2e (X'{t) - X''(0) + DpH\P)tj {I + DI^u'{x.')) dwt + e V + £'kw'(x^)| 
Hence, 



dt 



dt. 



E 



X^(T) - X^(0) + DpH {P)T 

T 



E 



^ E 



I 2e(x%t)-X%0)+DpH'{P)t){I + DJ,^u%x^})dwt+ f e^\I + D%^u%x^)\^ dt 

re'\I + Dl^u^^n\'dt 
Jo 



Dividing by T and letting T go to infinity 



lim E 



X'^(r) - X'^(O) + DpH {P)T 
T 



lim E 



T ^2 



T 



■dt 



Dj,^u''\''d0f,, 



2 / \I + Dl^u'\^d9^.<2ne^ + 2e^ ' "^ .,s,2 

T" Jt 

2 I o / l7-i_„.e|2ja I o / I n IT 7-i„Tj'^|2 



<2ne^ + 2 \Dpu^-' dO^. +2 \DpH - DpH \^ dO^., 

where we used (4.2). 
D 

We conclude the section with a proposition which shows in a formal way how 
much relation (1.3) is "far" from being an actual change of variables. Let us set 
'w^{x,P) :~ P ■ X + u^{x,P), where u^{x,P) is a ^"-periodic viscosity solution of 
(1.17), and let A; € Z". We recall that in the convex setting the following weak 
version of the change of variables (1.3) holds [EGOl, Theorem 9.1]: 

lim/ ^{d'},u{x,P)) d0f,^ ^{X)dX, 

for each continuous Z"-periodic function $ : M" -^ M, where 

' u{x, P + hei) — u{x, P) u{x, P + hcn) — u{x, P) 



D%u{x,P) := 



h 



h 
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ei, . . . , e„ being the vectors of the canonical basis in K". The quoted resuh was proven 
by the authors by considering the Fourier series of $, and then analyzing the integral 
on the left-hand side mode by mode. The next proposition shows what happens for 
a fixed mode in the non convex case. 

Proposition 7.3. The following inequality holds: 



[k ■ DpH ) I e 



2-Kik-Dpw^ 



d9„ 



<27r|/c|2 (e^+ I \Dpu'fd0, 



^- , , \DpH-DpH \Ue^. 



Proof. Recalling identity (3.8) with 



^(x) - e^'^^'^-^^'^^f"^^) 



we obtain 
= 



e.P„2TTik-Dpw^ 



L^^^e 



de„ 



2'Ki I e 

IT 

2'Ki I e 



2-!Tik-Dpw^ [re.P 



[L^'^ {k ■ Dpw^) - Trie^\D^{k ■ Dpw^)\^] dOf, 



27Tik-Dpw^ 



k ■ DpH' - Trie^\D^{k ■ Dpw 



e\\2 



d6,. 



where we used (4.9) and the fact that w^ — P ■ x + u^. Thus, thanks to estimate (4.2) 



(k ■ DpH 



<27r|/c|^ 



< 27r|/c|^ 



JIirik-Dpw^ 



de., 



<ne' \D,{k-Dpw')\'de^. 



\Di,u^U9^. 



\Dpu'\''d9^, 



\DpH-DpH Yde^,. 



D 

Remark 7.4. When H is uniformly convex, thanks to (4.5) the last chain of 
inequalities becomes 



[k-DpH) I e 



2-Kik-Dpw^ 



d9„ 



< C\k\^e^ (1 + trace {DlpH')) . 



Thus, if trace {DppH ) < C, the right-hand side vanishes in the limit as e — > 0, and 
we recover [EGOl, Theorem 9.1]. 

8. Compensated compactness. In this section, some analogs of compensated 
compactness and Div-Curl lemma introduced by Murat and Tartar in the context of 
conservation laws (see [Eva90], [Tar79]) will be studied, in order to better understand 
the support of the Mather measure /i. Similar analogs are also considered in [EvalO], 
to investigate the shock nature of non-convex Hamilton- Jacobi equations. 
What we are doing here is quite different from the original Murat and Tartar work (see 
[Tar79]), since we work on the support of the measure 0^e. Besides, our methods work 
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on arbitrary dimensional space R" while usual compensated compactness and Div- 
Curl lemma in the context of conservation laws can only deal with the case n — 1,2. 
However, we can only derive one single relation and this is not enough to characterize 
the support of n as in the convex case. To avoid confusion, when necessary we will 
explicitly write the dependence on the P variable. 

Let be a smooth function from T" x M" ^ R, and let p^ = {(/), H}0f,. + 

-7r4>pjPk^%ix-^xiXk^t^'- By (3.7) and (4.1), there exists C > such that 

/ \p'\dx<C. 

So, up to passing to some subsequence, if necessary, we may assume that p^ ^ p as 

a (signed) measure. 

By (5.2), /9(T") = 0. We have the following theorem. 

Theorem 8.1. The following properties are satisfied: 
(i) for every (j) € C(T" x R") 

DpH -{p- P)(l){x,p)dn^ udp; (8.1) 

(a) for every (f) E C(T" x R") and for every rj E C^{T'^), 

DpH ■ Dr)(j){x,p)dp^ rjdp. (8.2) 

T"xR" JT" 

Proof. Let w'^ — (f>{x, P + D^u^). Notice first that 
/ DpH -{p- P) (j}{x, p) dp = lim / DpH{x, P + Dxu") ■ OxU^w^dO^. . 

JT"xR" '^-^"JT" 

Integrating by parts the right hand side of the above equality we obtain 
/ DpH{x, P + Dxu") ■ OxU^w^dep. = - / u''div{DpHw''e^.)dx 

= - / u%div{DpH9p.)w' + DpH ■ Dxw'e^.)dx:^ I u'{ — Aep.w''-DpH-Dxw'0^ 
After several computations, by using (2.1) we get 



)dx. 



DpH ■ Dxw' = -{(/), H} + -^pA< 



Hence 



2 2 2 

^Ae^.w' - DpH ■ Dxw^e^,. ^ ^—M^.w' + {(i>,H}e^. - ^—4>p^Au%fi^,s 

2 2 2 

^yAw^e^. + y(div(i^,0^.zi;^) - d:w(Dxw'e^.)) + {0,ff}0^. - ^—4>pAu%^ep 



„ '-i'^r. .^ . f-l'rf 



+ 4>p,xXx,X, + 'i>X^X, + (jjpAKj^t^ 



n \^PiPk"-XiXj"-XiXk ' 'rpjXi^XjXi I 'TXiXi i Ypi^^XiJ")^ 

-+- ^^^ . ^^^ .J. ^ p- ^ ^ ^^ . ^^^ J. ^ ^ p- J J I ^^.^,^^j^^^ 

£2 £2 £2 

-^4>x,x,d^,.+ — (l)p^x,U%^^^ + - 



2 2 

+ ^-{AiY{DJp.w^') - dW{Dxw'ep.)) + {^,H}e^. - ^0p^A<^V 

£2 £2 

=9" + —4>xix,d^,s +—(j)p^^^u%.^^ + —{A\v{Dx9p.w^) - div{DxW^6p.)). 
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Therefore 

DpH-{p-P)c^{x,p)dn 



T"xl 



lim 

e->0 



P +^ 



„ YXiXi^fJ.'^ 2 ^Pi^i XjXi 



— {div{Dxe^.w^)-div{D^w^0f,s)) dx. 

(8.3) 



Since u^ converges uniformly to u, 



lim 



u^p^dx — udp 



The second term in the right hand side of (8.3) obviously converges to as £ — )■ 0. 
The third term also tends to by (4.1). 
Let's look at the last term. We have 



lim — 

lim — / dYv{Dxu'^w^)0^ 

£^0 2 Jt" 



lim 



>o 2 jj 

-2 



-D^u'^ ■ Dj.Onew'^ + Dr^u"^ ■ D^w^Onedx 



e^ r 
lim — / {Au^w^ + 2D^u'' ■ D^w" 
E^o 2 Jt" 



u'^{Aw{Dri.0^cw^)—diY{Dxw'^6^E))dx 
D^u'^ ■ DxW^Ofj^cdx 
< lim Ce^ [ iDl^u'W^.dx < lim Ce = 0, 

which implies (8.1). Relation (8.2) can be derived similarly. D As a consequence, we 
have the following corollary. 

Corollary 8.2. Let m(-,P) be a classical solution of (1.2), and let p, he the 
corresponding Mather measure given by Theorem 1.2. Then, 

DpH ■ [p — P — Dxu) = in supp p. 

Proof. By (8.1) and (8.2) 



/ DpH -{p-P - Dxu)(t>diJL^O, 



for all (f). Therefore, the conclusion follows. D 

9. Examples. In this section, we study non-trivial examples where the Mather 
measure ji is invariant under the Hamiltonian dynamics. Notice that, by (5.2), the 
Mather measure /i is invariant under the Hamiltonian dynamics if and only if the 
dissipation measures (jnkj) vanish. An example in Section 10 shows that this is not 
always guaranteed. As explained in [EvalO], the dissipation measures nikj record the 
jump of the gradient D^u along the shock lines. 

We investigate now under which conditions we still have the invariance property 
(1). We provide some partial answers by studying several examples, which include 
the important class of strongly quasiconvex Hamiltonians (see [FS05] ) . 

9.1. H is uniformly convex. There exists a > so that Dp H > a > 0. 
Let A = in (5.5) then 

= / -ffpfcPj dm^j , 

JT"xR" 

which implies rukj = for all 1 < k, j < n. We then can follow the same steps as in 
[EGOl] to get that p also satisfies (2). 
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9.2. Uniformly convex conservation law. Suppose that there exists F{p,x), 
strictly convex in p, such that {F, H}^0. Then m = 0. 

9.3. Some special non-convex cases. The cases we consider here are somehow 
variants of the uniformly convex case. 

Suppose there exists (/) uniformly convex and a smooth real function / such that 
either (p = .f(H) or H = f{<i>). Then, by (5.2) we have mkj = for all fc, j. In 
particular, if _ff = /(</)) with / increasing, then H is quasiconvex. 
One explicit example of the above variants is H{x,p) = (|pp + V(x))'^, where V : 
T" — > M is smooth and may take negative values. Then H{x,p) is not convex in p 
anymore. Anyway, we can choose (j){x,p) — |pp + V{x), so that H{x,p) — (0(x,p))^ 
and (j) is uniformly convex in p. Therefore, /x is invariant under the Hamiltonian dy- 
namics. 

9.4. The case when n — 1. Let's consider the case H{x,p) — H{p) + V{x). 
In this particular case, property (H3) implies that \H{x,p)\ — J- oo as \p\ — > +oo. Let 
us suppose that 

lim H{p) = +CX). 

|p|^+oo 

Assume also that there exists po G K such that H'{p) = if and only if p = Po 
and H"{po) ^ 0. Notice that -ff(p) does not need to be convex. Obviously, uniform 
convexity of H implies this condition. 

We will show that m\\ = 0, which implies that /i is invariant under the Hamiltonian 
dynamics. From our assumptions, we have that H'{p) > for p > po, H'{p) < for 
p < Po and hence H"{po) > 0. Then there exists a neighborhood (po — r, po + ^) of po 
such that 

H"{p)>^^^, ypeip,-r,po + r). 

And since the support of mn is bounded, we may assume 

supp(mii) C Tx [-M,M], 

for some M > large enough. We can choose M large so that {po — r,po + r) C 

(-M,M). 

Since |i?'(p)P > for p e [-M, M] \ {po -r,po + r) and [-M, M] \ {po -r,po + r) is 

compact, there exists 7 > such that 

|ff'(p)|'>7>0, Vpe [-M,M]\(po-r,po + r). 

Hence, by choosing A ;» 

X\H'{p)\' + H"{p)>^^^, Vpe[-M,M], 

which shows mn = by (5.5). 

9.5. Case in which there are more conserved quantities. Let's consider 

H{x,p) = H{p) + V{xi + ... + Xn), 
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where F : T ^ M is smooth. 

For k 7^ j, define $'^-' = Pk — Pj- It is easy to see that {H, ^'^^j = for any k 7^ j. 

Therefore {if, ($*'^)^} = for any k ^ j. 

For fixed k ^ j, let (j) = ($'=^)2 in (5.2) then 



T"xl 



("T-fefc ~ 2TOfej + nijj) dxdp = 0. 



The matrix of dissipation measures {mkj) is non-negative definite, therefore mkk — 
2mkj + rrijj > 0. Thus, mkk — '^rnkj + rnjj = for any k ^ j. 

Let e G (0, 1) and take £_ == (^1, ...,^„), where £_k = 1 + £jCj = ^1 a-^d ^i = otherwise. 
We have 

< TUkj^k^j = (1 + e)'^mkk - 2(1 + e)mkj + rrijj = 2e{mkk - mtj) + e^mkk- 

Dividing both sides of the inequality above by e and letting e — > 0, 

mkk - mkj > 0. 

Similarly, tujj — mkj > 0. Thus, mkk — m,kj = rrijj — rrikj = for all k ^ j. 
Hence, there exists a non-negative measure m such that 

mkj ~ m>0, V k,j. 

Therefore, (5.5) becomes 

We here point out two cases which guarantee that tti = 0. In the first case, assuming 
additionally that H{p) ~ Hi(pi) + ..._ff„(p„) and H2,-.-,Hn are convex, but not 
necessarily uniformly convex (their graphs may have flat regions) and Hi is uniformly 
convex, then we still have m — 0. 

In the second case, suppose that H{p) = H{\p\), where H : [0, 00) — > R is smooth, 
iJ'(O) = 0,7J"(0) > and H'{s) > for s > 0. Notice that H is not necessarily 
convex. This example is similar to the example above when n — 1. Then for p 7^ 

and at p = 

A(^i/p^(0))'+^i/p^p,(0) -ni7"(0) >0. 

So, we can choose r > 0, small enough, so that for \p\ < r 

j 3,k 

Since the support of m is bounded, there exists M > large enough 

suppmcT" X {p: \p\ < M}. 
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Since min^gj^ m] H'{s) > 0, by choosing A ^ 0, wc finally have for \p\ < M 

f n ri . /„„,„. niin^g[ M]^'(s) 
for /5 = — mm -^ ii (0), 



2 I M 

Thus m — 0, and therefore /i is invariant under the Hamiltonian dynamics. 

9.6. Quasiconvex Hamiltonians: a special case. Let's consider 

H{x,p)^H{\p\) + V{x), 

where i? : [0, oo) ^ M is smooth, i7'(0) = 0, iJ"(0) > and H'{s) > for s > 0. 
Once again, notice that H is not necessarily convex. We here will show that {nrijk) ~ 0. 
For p 7^ then 

{\H,^H,, + H,^,,)m,u = ^(mn + ... + m„„) + (^\(H'f + iJ" - ^^ ^^^. 

For any symmetric, non- negative definite matrix m = (fnjk) we have the following 
inequality 

< PjPkrrijk < Ipptracem = |pp(mii + ... + m„„). 

There exists r > small enough so that for |p| < r 



¥\ > >^'^-' 






< \h"{0). 



Hence for \p\ < r 

(XHp^Hp^ + Hp^pjnijk > -H"{0){mii + ... + m„„). 

Since the support of {rrijk) is bounded, there exists M > large enough 
supp mjk C T" X {p : IpI < M}, V j, k. 

Since min^^gr^ ^i H'{s) > 0, by choosing A 3> we finally have for \p\ < M 
{XHp.Hp^ + Hp^p^)Tnjk > l3{mii + ... + to„„), 

for p — mm < , — — ■ > > 0. 

^ 12' M J 

We then must have rtiw + ...to„„ — 0, which implies (jrijk) — 0. Thus, /i is invariant 

under the Hamiltonian dynamics in this case. 

We now derive the property (2) of ^ rigorously. Since the support of ji is also bounded, 

we can use a similar procedure as above to show that <f){x,p) = e^^^^'^^ is uniformly 

convex in T" x i?(0, M) D supp(/i) for some A large enough. 

More precisely, 

0p,P.Oa > e^^m', e e M", {x,p) e T" X B(0,M), 
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for /? chosen as above. Then doing the same steps as in [EGOl], we get /i satisfies (2). 
There is another simple approach to prove (2) by using the properties we get in this 
non-convex setting. Let's just assume that u is C^ on the support of /i. 
By Remark 8.2, it foUows that DpH.{p — P — Du) = on support of ji. And since 

DpH{x,p) = H'{\p\)y- for p 7^ and H'{\p\) > 0, we then have p.{p - P - Du) = 

on support of /i. Hence \p\'^ ^ p.{P + Du) on supp(/i). 

Besides, H{x,p) — H{x,P + Du{x)) = H{P) on supp(/i) by property (a) of Mather 
measure and the assumption that u is C^ on supp(/i). It follows that iJ(|p|) = 
H{\P + Du\). Therefore, \p\ = \P + Du\ by the fact that H{s) is strictly increasing. 
So we have |pp = p.(P+Du) and \p\ — \P+Du\ on supp(/x), which implies p — P+Du 
on supp(/j.), which is the property (2) of fi. 

9.7. Quasiconvex Hamiltonians. We treat now the general case of uniformly 
quasiconvcx Hamiltonians. We start with a definition. 

Definition 9.1. A smooth set A C M" is said to be strongly convex with convexity 
constant c if there exists a positive constant c with the following property. For every 
p G dA there exists an orthogonal coordinate system (qi, . . . , Qn) centered at p, and a 
coordinate rectangle R = (ai, 6i)x . . . x (a„, 6„) containing p such that TpdA ~ {qn = 
0} and A n i? C {(J e i? : c^"J^^ \qi\^ < Qn < &n}- The previous definition can be 
stated in the following equivalent way, by requiring that for every p G dA 

(Bpv) • V > c|vp for every v e TpdA, 

where Bp : TpdA x TpdA ^ R is the second fundamental form of dA at p. 

We consider in this subsection strongly quasiconvex Hamiltonians. That is, we 
assume that there exists c > such that 

(j) {p G T" : H{x,p) < a} is strongly convex with convexity constant c for every 
a € R and for every a; € T". 
In addition, we suppose that there exists a e M such that for every x E T" 
(jj) There exists unique p e M" s.t. DpH{x,p) — 0, and 

DlpH{x,p) > a. 

Notice that the special case presented in Section 9.6, where the level sets are spheres, 
fits into this definition. We will show that under hypotheses (j)-(jj) there exists A > 
such that 

A DpH €) DpH + DpH is positive definite. 

From this, thanks to relation (5.5), we conclude that nikj ~ 0. First, we state a 
well-known result. We give the proof below, for the convenience of the reader. 

Proposition 9.2. Let (j)-(jj) be satisfied, and let ix*,p*) e T" x M" be such 
that DpH{x*,p*) ^ 0. Then 

DpH{x*,p*)±Tp,C and DlpH{x* ,p*) = \DpH{x* ,p*)\Bp, , (9.1) 

where Bp* denotes the second fundamental form of the level set 

C := {p e R" : H{x*,p) = H{x*,p*)} 

at the point p* . 
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Proof. By the smoothness of H, there exists a neighborhood U C K" of p* and 
n smooth functions v : U ^^ 5"^^, Ti : U ^^ 5"^^, z = 1, . . . ,n — 1, such that for 
every p G U the vectors {ti{p), . . . , T„_i(p), i^(p)} are a smooth orthonormal basis of 
K", and for every p € U (IC Ti{p), . . . ,r„_i(p) G TpC. Let now i, j G {1, . . . , n — 1} 
be fixed. Since 

H{x*,p)^a ypeU, 

differentiating w.r.t Ti{p) we have 

DpH{x*,p)-n{p)^0 VpGC/nC. (9.2) 

Computing last relation ai p = p* we get that DpH{x* ,p*) _L Tp*C. Differentiating 
(9.2) along the direction Tj{p) and computing at p = p* 

{DlpH{x*,p*)T,{p*)) ■ n{p*) + DpH{x*,p*) ■ (Dpn{p*)T,{p*)) = 0. (9.3) 

Notice that by differentiating along the direction Tj{p) the identity Ti{p) ■ v{p) = 
and computing at p* we get 

{Dpn{p*)Tj{p*)) ■ v{p*) = - {Dpv{p*)Tj{p*)) ■ n{p*). 

Plugging last relation into (9.3), and choosing v{p*) oriented in the direction of 
DpH{x*,p*) we have 

{DlpH{x*,p*)T,{p*)) ■ n{p*) ^ -\DpH{x*,p*)\ {Dpn{p*)T,{p*)) ■ i^ip*) 

= \DpHix*,p*)\ {Dpv{p*)T,{p*)) ■ nip*) ^ \DpH{x\p*)\ (Bp.T,(p*)) • n{p*). 

D 

For every vector v G M", we consider the decomposition 

V = i;||v" + f_LV , 

with ull,w^ G R, |vll| = |v-L| = 1, v" G TpX, and v^ G [Tp.C)^ . By hypothesis (jj) 
and by the smoothness of iJ, there exist r > and a' G (0,a), independent of {x,p), 
such that 

DlpH{x,p) > a' for every {x,p) G {\DpH\ < t}. 

Let us now consider two subcases: 
Case 1: {x,p) G {\DpH\ < t} 

First of all, notice that 

XDpH ® DpHv ■ V = X\DpH ■ v\^ = \v\ \DpH\^ . 
Then, we have 

(A DpH ® DpH + DlpH)v ■ v = \v\ \DpHf + {DlpHv ■ v) > a>|2. 
Case 2: {x,p) G {\DpH\ > r} 
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In this case we have 



DIj,Hw\\ ■ v" > c\DpH\ 



which then yields 



> cvfADpH\ + 2v\\v^{Dl Hw^\ ■ v^) + vI{DI Hv^ ■ v-^) 



By (5.3) we have 



Dlr,H\ < C along supp/x. 



Thus, 



(A DpH ® DpH + DlpH)v ■ v 

>\vl \DpH\^ + cvl\DpH\ + 2w||i;_L(£'ppi?vll • v-^) + vl{DlpHw^ ■ v^) 
> v'i (A \DpH\'^ -C)- 2C\v\\\\v^\ +cvl\DpH\ 

>wi(Ar2-c(l + ^)) +vl{cT-CTf). 



Choosing first rf < ^r, and then 



A > — 1 + — 



we obtain 



(A DpH ® DpH + Dl H)v ■ v > a"\v\^, 



for some a" > 0, independent of {x,p). 
General Case 

In the general case, we have 

{XDpH ® DpH + DlpH)v ■ v > "f\v\^, 

where 7 := min{a', a"}. 

Similar to the case above, we basically have that (p{x,p) — e^^^^'^' is uniformly convex 
on the support of fi for A large enough. Hence, by repeating again the same steps as in 
[EGOl], we finally get that /i satisfies (2). As already mentioned in the introduction, 
we observe that one could also study the case of uniformly convex Hamiltonians by 
duality, that is, by considering a function $ : M ^ K such that $(_ff (x, •)) is convex 
for each x G T". In this way, the dynamics can be seen as a reparametrization of the 
dynamics associated to the convex Hamiltonian ^{H). 

10. A one dimensional example of nonvanishing dissipation measure to. 

In this section we sketch a one dimensional example in which the dissipation measure 
TO. does not vanish. We assume that the zero level set of the Hamiltonian iJ : T x M — >■ 
R is the smooth curve in Figure 10.1, and that everywhere else in the plane {x,p) the 
signs of H are as shown in the picture. In addition, H can be constructed in such a 
way that {D^H,DpH) ^ (0,0) for every {x,p) e {{x,p) e T x M : H{x,p) ^ 0}. That 
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H{x,p)>Q : 


— 9{x) 
-{H{x 



Fig. 10.1. {H{x,p) = 0}. 



Fig. 10.2. g{x). 



is, the zero level set of H does not contain any equilibrium point. Consider now the 
piecewise continuous function g : [0,1] — > M, with g(0) ~ .9(1), as shown in Figure 
10.2. Then, set 



P := I g{x) dx, 
'o 



u{x, P) := -Px + I g{y) dy 



and define 



One can see that m(-, P) is the unique periodic viscosity solution of 

H{x,P + D,u{x,P))^0, 

that is equation (1.2) with H{P) = 0. Assume now that a Mather measure /i exists, 
satisfying property (1). Then, the support of fi has necessarily to be concentrated 
on the graph of g, and not on the whole level set {H = 0}. However, any invariant 
measure by the Hamiltonian flow will be supported on the whole set {H = 0}, due 
to the non existence of equilibria and to the one-dimensional nature of the problem, 
thus giving a contradiction. 
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Rezakhanlou for very useful discussions on the subject of the paper. 
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